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=3 ■ Abstract 



We propose a practical entanglement classification scheme for general multipartite 
pure states in arbitrary dimensions under local unitary equivalence by exploiting the 
^ i high order singular value decomposition technique and local symmetries of the states. 

By virtue of this scheme, the method of determining the local unitary equivalence of 
n-qubit states proposed by Kraus is extended to the case for arbitrary dimensional 
multipartite states. 

in 

O ■ PACS numbers: 03.67.Mn, 02.10.Xm, 03.65.Ud 



^ I Entanglement is one of the most extraordinary features of quantum theory. It lies at 

■ the heart of quantum information theory and is now regarded as a key physical resource in 
realizing many quantum information tasks, such as teleportation and quantum computation, 
etc. pp. In practice, people may confront various forms of entangled states even in one single 
physical system. Though superficially showing up with different features, usually not all these 
entangled states are functionally independent; they may be intrinsically the same where the 
entanglement property is concerned. Two entangled states are said to be equivalent in 
implementing the same quantum information task if they can be obtained with certainty 
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from each other via local operation and classical communication (LOCC). Theoretically, 
this LOCC equivalent class is defined such that within this class any two quantum states 
are interconvertible by local unitary (LU) operators [2J. If we release the constraint of the 
local unitary operator to invertible local transformation, the LOCC equivalence then turns 
to the widely discussed stochastic local operation and classical communication (SLOCC) 
equivalence. 

The classification of a bipartite system under LU can be done by singular value (Schmidt) 
decomposition. For the classification of a pure multipartite entangled state, a canonical 
method was proposed in Ref. [3], though it was given only in a set of constraints on the 
coefficients of quantum state. Later, this method was reformulated into a compact form 
where the existence of SLOCC normal forms was also proved. Recently, by introducing the 
standard form for multipartite states, Kraus proposed a general way to determine the local 
unitary transformation between two LU equivalent n-qubit states [5]. Even though these 
beautiful results have been obtained, a generic method or methodology is still missing in the 
classification of general multipartite entangled states [6]. 

Much research indicates that symmetry study is very helpful in the classification of a 
general multipartite system. It is broadly realized that making restrictions on the quantum 
states, i.e., satisfying some symmetric properties, is a feasible way in the entanglement 
classification [TJIH]. By virtue of the Majorana representation, a lot of effort was spend on the 
study of the relation between permutation symmetry and the classification of a multipartite 
entangled state under both LOCC and SLOCC [9l[T0l[TT]. It was recently realized that there 
exists one kind of local (internal) symmetry under SLOCC equivalence for a continuously 
entangled state [12]. Unlike the permutation symmetry among different partites, the internal 
symmetry is of a localized one. In Refs. [131 [H US] , people investigated the nature of the 
n-qubit entangled state via analyzing the stabilizers of the symmetric state, which can also 
be regarded ClS db local unitary transformation. 

In this work, we present a new and easy to operate entanglement classification method 
under the LU equivalence, by virtue of the technique of high order singular decomposition 
(HOSVD) [ini [13 and by exploiting the local symmetries of the state [12] . The procedure 
of this method goes as follows: Express the quantum state into a complex tensor, analyze 
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the possible local symmetry property of the quantum state, apply the HOSVD on the tensor 
form of the quantum state, and then fully classify a general multipartite state into LU 
inequivalent classes by sequential use of HOSVD and local unitary symmetry. 

A general TV-partite entangled quantum state in dimension /i x J2 x • • • x /jv can be 
formulated as the following form: 

Il,l2,-,lN 

I*) = V'ni2...i;vl^l)N2)-|«iv) , (1) 

where £ C are coefficients of the quantum state in representative bases. Two 

quantum states are said to be equivalent if they are interconvertible by a certain type of 
operators, which can be schematically expressed as 

N 

1^') = (g)t/«|^) 
i 

i' i' ■■■ i' 

<i2-ijvKl'^2,--- ,^iv) • (2) 



«1,«2,--- ,IN 



Here, the coefficients ipij^i2...iM can also be treated as the entries of a tensor and hence the 
quantum states can be represented by high dimensional complex tensors. From the tensor 
form of ^, the operator t/*^") acting on the nth partite is defined as 

([/*• ^^)iii2— in-l«J,2n+l---jjv ~ ^»l«2---in-12nin+l— jjv'^i^]„ • (^) 

In case J/^'^^'s are unitary operators, quantum state ^' is said to be LU equivalent to ^, i.e., 

^' = 0;7W^. (4) 

i 

Suppose a concerned quantum state is inherited with local symmetries, like 

ipW^ = ^, (5) 
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and is always assisted with classical communication hereafter in our discussion; the local 
symmetry group G = {g\g = (S^n-^''"'*} then forms the stabilizer of \1'. Here P^"-* is the 
operator that acts on the nth partite. From Eq.(jl]), we can get 

^Wp«^Wt^' (6) 



with [/^*^ the unitary operators. Clearly, every local unitary equivalent state has the local 
symmetry of P^'''')'{^k U^^^^), which is isomorphic to that of ^. Consequently, 

we have the following proposition. 

Proposition 1 // \E'' is local unitary equivalent to \E', then the stabilizers of the local sym- 
metries of the two states are unitarily equivalent. 

By taking = a|ll) + ^|22) + 7|33) +7|44), where |ap + |/3|2 + 2|7|2 = 1 as an example, 
this 4x4 entangled state can be expressed in a matrix form: 



4x4 



/a 0\ 

0/300 

7 

\0 7/ 



(7) 



It is invariant under the local transformation U^^^ = {e*^^ , e*^^ , u} and U^"^^ = {e , e , m*}, 
i.e., 

^(l).^.(t/(2))T = ^_ (8) 

The ^^4x4 is then a f/*^^-* ® ^7*-^^ invariant quantum state. 

To study the LU classification of a multipartite system, it is convenient for us to use the 
notation of Ref.piiJ and to define the matrix unfolding of the tensor \E' G C^^^^ "^^ with the 
nth index as 
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Here is an J„ x (/„+iJ„_|_2 ■ ■ ■ iNhh ■ ■ ■ In-i) matrix. Considering the high order singular 
value decomposition developed in Ref. |16], there exists a core tensor S for each \E', that is, 



^ = f/(i) ® [/(2)...®[/Ws , (10) 

where S forms a same order tensor in the Hilbert space Ii x I2 x ■ ■ ■ x Ij^. Any order 
tensor Si^=i, obtained by fixing the nth index to i, has the following properties [16] : 

(Si„=i, S,„=j) = 6ij ^ ^ , (11) 

with cTj^""* > crj"^ and \/ i < j for all possible values of n. Here, the singular value 
symbolizes the Frobenius norm cr-"^ = ||Sj^=j|| = ^J T,i^=i), where the inner product 



From the definition, we know that if the quantum state is LU equivalent to \E', then they 
can be transformed into the same core tensor S . However, because there is no information 
on whether different core tensors are local unitary equivalent or not, S has nothing to do 
with the entanglement classification of the multipartite state yet. Hence, to achieve the 
entanglement classification from this point, the following analysis is necessary. 

Since every unitarily equivalent quantum state possesses the same local symmetric prop- 
erty, we can pick up a typical state to represent one equivalent class. A natural choice is 
the core tensor for each class. In the matrix unfolding form of the core tensor, S(„) can be 
expressed as 

^(n) = [/("^ ■ ■ ([/("+^) ® ■ ■ • ® U^^^ ■ [/(^) ■ ■ ■ [/("~i))T , (12) 

where f/*^"^ is constructed by the left singular vectors of "^(n)- 

There exist such quantum states \E' that some of the singular values of their unfolding 
matrix \1'(„) are identical, which is similar to the singular value 7 in the matrix case in Eq.(I7]) 
and is also stated in Property 4 of Ref. |16j. The \E'(„) has the following property: 

= [/WpW . • (^7("+^) ® ■ ■ • ® ?7(^) ® [/(^) ■ • ■ [/("-i))T , (13) 
where P^"-* is a block-diagonal matrix consisting of unitary blocks with the same partitions 
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as that of the identical singular values, and S'^^^ = P("^'I"S(„). For such states we have 

N 

S' = (g)p(")ts . (14) 

n=l 

Here, P'-"^ = diag{e*^^, ■ ■ ■ , e*^^"} in case there is no identical singular value for Thus, 
if S is a HOSVD of \E', then the S' is as well. In other words, the HOSVD of possesses 
local symmetry properties of <S)n=i -P*-"^^- 

From the above discussion, we may conclude that if we introduce the symmetry prop- 
erty to compensate the nonuniqueness of the core tensor, it will then serve as the unique 
representative state for each local unitary equivalent class. We call the representative state 
S associated with specific symmetry group <S)n=i ^^^^ unique canonical form of the en- 
tanglement class under LU equivalence. Two quantum states are LU equivalent if and only 
if they have the same core tensor (representative state) or their core tensors are related by 
the symmetry group <S>n=i P^""^- 

Proposition 2 The high order singular value decomposition of multipartite state is a LU 
classification of entanglement up to the local symmetry ^^=1 -P*'"''- 

Proof: First, it is easy to see that if two quantum state are LU equivalent, i.e., \E'' = ^f^U^^^"^, 
then they can be transformed into the same core tensor E. Second, suppose two core tensors 
are related by local unitary transformations S' = V^^'^T,, where V^""^ are different 
from ^n=i -P*-"''; then by taking the nth unfolding as an example, we have 

This means that the singular values are not uniquely determined, which is in contradiction 
with Property 4 of Ref. [16]. Therefore, we can conclude that the core tensor S up to the 
local symmetry ^^=1 ^^"^^ unique, or, in other words, the core tensor of the quantum 
state can be treated as the LU classification of a multipartite quantum state up to the local 
symmetry (g)^^iP("). QED. 

Now, for a general arbitrary dimensional A^-partite pure state, the canonical form can 
be constructed in the following way. First, express the tensor form quantum state in matrix 
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unfolding form [see Eq.([9])]. Then, compute the singular value decomposition of the 
matrix \E'(„) = U^"''^ AV^'^\ Here unitary matrices U^"'\ composed by left singular vectors 
of ^E'(n)5 for N indices will be obtained, and this is a simple process while only matrix 
operations are involved. Finally, the core tensor is constructed by 

N 

S = (g) [/(")t^ . (16) 

n=l 

The representative state S with its associated sjTiimetry group in the form of Eg. (1141) will 
serve as the unique canonical form of 

Given two canonical forms, i.e., the obtained representative states S and S with some 
local symmetries of Q*-"^-* and <^„ Q^"'\ respectively, the procedure of verifying the ex- 
istence of a local symmetry that connects the two representative states is straightforward. 
If singular values of these two representative states or the sizes of the unitary blocks of 
local symmetries are not the same, obviously there does not exist such a local symmetry 
that the two tensors can be transformed into each other. Otherwise, the existence of the 
local symmetry is equivalent to the existence of solutions to the equation S = (^^Q'-^^S 
(or S = (S^n^''"^^)' "whsre the elements in S and S are known quantities whereas they are 
unknown variables in Q^") By virtue of the mathematical technique relinearization 

[13, [IS], one can readily know whether the solutions to the equations exist or not (practical 
examples are given in Ref. [20]). 

Last, we show briefly the relation of the core tensor and Kraus's canonical form l^&st) in 
Ref. [5|. The sorted trace decomposition of the nth partite requires Tr^j,^(|\E'st)(\l'st|) to be 
diagonalized. For a general quantum state 

= 1^"-^ ^ili2---in-ijnjn+i---iiv Nl^2 ' ' ' ' ' ' ^Af) , (17) 

the requirement of trace diagonalization asks 

AS")>A5")>--->Ai")>0. (18) 
7 



Here, A}"^ = E^j <i2-i„-iji„+i-iiv«ni2-*n-iiWi-*iv is just the Frobenius norm ||(^st)i„=ilP 
in high order singular decomposition employed in our classification approach. 

In summary, by virtue of the symmetry property study and tensor decomposition we 
propose an effective and easy to operate method for the local unitary classification of a 
multipartite entangled state, which splits the mathematically difficult problem into different 
relatively easier ones. That is, we first play tensor decomposition on the quantum state and 
obtain the quasinormal form (core tensor), where many different quasinormal forms may 
actually belong to the same entanglement class; then, by exploiting the local symmetries 
between the quasinormal forms, we arrive at a uniquely defined entanglement class. By this 
approach, the verification of LU equivalence of n-qubit states proposed by Kraus is extended 
to the case for arbitrary dimensional multipartite states in a more easily understandable 
and computable algorithm. From a methodological point of view, the approach in this work 
is different from previous ones in the literature, where normally the entangled state is first 
divided into apparently inequivalent coarse-grained sets under certain criteria (like range [2T], 
partition number [8], ranks [22], etc.) and then more fine-grained classification procedures 
are performed. This new approach might also be useful in the study of classification under 
SLOCC or other characters of the general multipartite entangled state. Finally, it should be 
mentioned that the tensor decomposition method may also play an important role in other 
entanglement property studies [23] . 
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